A three-dimensional immersed boundary method for non-Newtonian fluids

Fluid-structure-interaction (FSI) phenomenon is common in science and engineering. The fluid involved in an FSI problem may be non-Newtonian such as blood. A popular framework for FSI problems is Peskin's immersed boundary (IB) method. However, most of the IB formulations are based on Newtonian fluids. In this letter, we report an extension of the IB framework to FSI involving Oldroyd-B and FENE-P fluids in three dimensions using the lattice Boltzmann approach. The new method is tested on two FSI model problems. Numerical experiments show that the method is conditionally stable and convergent with the first order of accuracy

An IB Method for Non-Newtonian-Fluid Flexible-Structure Interactions in Three-Dimensions

Problems involving fluid flexible-structure interactions (FFSI) are ubiquitous in engineering and sciences. Peskin’s immersed boundary (IB) method is the first framework for modeling and simulation of such problems. This paper addresses a three-dimensional extension of the IB framework for non-Newtonian fluids which include power-law fluid, Oldroyd-B fluid, and FENE-P fluid. The motion of the non-Newtonian fluids are modelled by the lattice Boltzmann equations (D3Q19 model). The differential constitutive equations of Oldroyd-B and FENE-P fluids are solved by the D3Q7 model. Numerical results indicate that the new method is first-order accurate and conditionally stable. To show the capability of the new method, it is tested on three FFSI toy problems: a power-law fluid past a flexible sheet fixed at its midline, a flexible sheet being flapped periodically at its midline in an Oldroyd-B fluid, and a flexible sheet being rotated at one edge in a FENE-P fluid

3D simulation of a viscous flow past a compliant model of arteriovenous-graft annastomosis

Hemodialysis is a common treatment for end-stage renal-disease patients to manage their renal failure while awaiting kidney transplant. Arteriovenous graft (AVG) is a major vascular access for hemodialysis but often fails due to the thrombosis near the vein-graft anastomosis. Almost all of the existing computational studies involving AVG assume that the vein and graft are rigid. As a first step to include vein/graft flexibility, we consider an ideal vein-AVG anastomosis model and apply the lattice Boltzmann-immersed boundary (LB-IB) framework for fluid-structure-interaction. The framework is extended to the case of non-uniform Lagrangian mesh for complex structure. After verification and validation of the numerical method and its implementation, many simulations are performed to simulate a viscous incompressible flow past the anastomosis model under pulsatile flow condition using various levels of vein elasticity. Our simulation results indicate that vein compliance may lessen flow disturbance and a more compliant vein experiences less wall shear stress (WSS)

A three dimensional implicit immersed boundary method with application

AbstractMost algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid. A drawback of such an explicit approach is a severe restriction on the time step size for maintaining numerical stability. An implicit immersed boundary method in two dimensions using the lattice Boltzmann approach has been proposed. This paper reports an extension of the method to three dimensions and its application to simulation of a massive flexible sheet interacting with an incompressible viscous flow

Two-dimensional Finite Element Model of Breast Cancer Cell Motion Through a Microfluidic Channel

A two-dimensional model for red blood cell motion is adapted to consider the dynamics of breast cancer cells in a microfluidic channel. Adjusting parameters to make the membrane stiffer, as is the case with breast cancer cells compared with red blood cells, allows the model to produce reasonable estimates of breast cancer cell trajectories through the channel. In addition, the model produces estimates of quantities not as easily obtained from experiment such as velocity and stress field information throughout the fluid and on the cell membrane. This includes locations of maximum stress along the membrane wall. A sensitivity analysis shows that the model is capable of producing useful insights into various systems involving breast cancer cells. Current results suggest that dynamics taking place when cells are near other objects are most sensitive to membrane and cytoplasm elasticity, dynamics taking place when cells are not near other objects are most sensitive to cytoplasm viscosity, and dynamics are significantly affected by low membrane bending elasticity. These results suggest that continued calibration and application of this model can yield useful predictions in other similar systems

Interactive 3D simulation for fluid–structure interactions using dual coupled GPUs

The scope of this work involves the integration of high-speed parallel computation with interactive, 3D visualization of the lattice-Boltzmann-based immersed boundary method for fluid–structure interaction. An NVIDIA Tesla K40c is used for the computations, while an NVIDIA Quadro K5000 is used for 3D vector field visualization. The simulation can be paused at any time step so that the vector field can be explored. The density and placement of streamlines and glyphs are adjustable by the user, while panning and zooming is controlled by the mouse. The simulation can then be resumed. Unlike most scientific applications in computational fluid dynamics where visualization is performed after the computations, our software allows for real-time visualizations of the flow fields while the computations take place. To the best of our knowledge, such a tool on GPUs for FSI does not exist. Our software can facilitate debugging, enable observation of detailed local fields of flow and deformation while computing, and expedite identification of ‘correct’ parameter combinations in parametric studies for new phenomenon. Therefore, our software is expected to shorten the ‘time to solution’ process and expedite the scientific discoveries via scientific computing

Modeling and Implementation of an Asynchronous Approach to Integrating HPC and Big Data Analysis

With the emergence of exascale computing and big data analytics, many important scientific applications require the integration of computationally intensive modeling and simulation with data-intensive analysis to accelerate scientific discovery. In this paper, we create an analytical model to steer the optimization of the end-to-end time-to-solution for the integrated computation and data analysis. We also design and develop an intelligent data broker to efficiently intertwine the computation stage and the analysis stage to practically achieve the optimal time-to-solution predicted by the analytical model. We perform experiments on both synthetic applications and real-world computational fluid dynamics (CFD) applications. The experiments show that the analytic model exhibits an average relative error of less than 10%, and the application performance can be improved by up to 131% for the synthetic programs and by up to 78% for the real-world CFD application

Simulation of blood flow past a distal arteriovenous-graft anastomosis at low Reynolds numbers

Patients with end-stage renal disease are usually treated by hemodialysis while waiting for a kidney transplant. A common device for vascular access is an arteriovenous graft (AVG). However, AVG failure induced by thrombosis has been plaguing dialysis practice for decades. Current studies indicate that the thrombosis is caused by intimal hyperplasia, which is triggered by the abnormal flows and forces [e.g., wall shear stress (WSS)] in the vein after AVG implant. Due to the high level of complexity, in almost all of the existing works of modeling and simulation of the blood-flow vessel-AVG system, the graft and blood vessel are assumed to be rigid and immobile. Very recently, we have found that the compliance of graft and vein can reduce flow disturbances and lower WSS [Z. Bai and L. Zhu, “Three-dimensional simulation of a viscous flow past a compliant model of arteriovenous-graft anastomosis,” Comput. Fluids 181, 403–415 (2019)]. In this paper, we apply the compliant model to investigate possible effects of several dimensionless parameters (AVG graft-vein diameter ratio Rgv, AVG attaching angle θ, flow Reynolds numbers Re, and native vein speed Vv) on the flow and force fields near the distal AVG anastomosis at low Reynolds numbers (up to several hundreds). Our computational results indicate that the influences of the parameters Rgv, θ, and Re lie largely on the graft and the influence of Vv lies largely on the vein. In any case, the WSS, wall shear stress gradient, and wall normal stress gradient and their averaged values on the graft are significantly greater than those on the vein. ACKNOWLEDGMENT

A deformable plate interacting with a non-Newtonian fluid in three dimensions

We consider a deformable plate interacting with a non-Newtonian fluid flow in three dimensions as a simple model problem for fluid-structure-interaction phenomena in life sciences (e.g., red blood cell interacting with blood flow). A power-law function is used for the constitutive equation of the non-Newtonian fluid. The lattice Boltzmann equation (the D3Q19 model) is used for modeling the fluid flow. The immersed boundary (IB) method is used for modeling the flexible plate and handling the fluid-plate interaction. The plate drag and its scaling are studied; the influences of three dimensionless parameters (power-law exponent, bending modulus, and generalized Reynolds number) are investigated

Scalings of Inverse Energy Transfer and Energy Decay in 3-D Decaying Isotropic Turbulence with Non-rotating or Rotating Frame of Reference

Energy development of decaying isotropic turbulence in a 3-D periodic cube with non-rotating or rotating frames of reference is studied through direct numerical simulation using GPU accelerated lattice Boltzmann method. The initial turbulence is isotropic, generated in spectral space with prescribed energy spectrum E(κ)~κm in a range between κmin and κmax. The Taylor microscale Reynolds number Reλ and Rossby number Ro are introduced to characterize the inertial, viscous, and rotational attributes of the system. The focus of this study is on the scalings of early inverse energy transfer and late energy decay in the development of turbulent energy under various conditions through combinations of m, κmin, κmax, Reλ and Ro. First, we demonstrate the validity of the simulation by confirming the quantitative dependence of the decay exponent n on the initial energy spectrum exponent m, at Reλ =255 and Ro=∞, varying the values of m, κmin and κmax. Second, at relatively low Reλ, the decay exponent for different initial spectra statistically fall in respective ranges, all of which agree well with the corresponding analytical predictions. Third, we quantitatively investigate the 3-D inverse energy transfer. Our findings include (i) the exponent of inverse energy transfer spectrum E(κ)~κσ depends on the initial spectrum exponent E(κ) ~ κm: if m<4, σ=m while if m≥4, σ=4; (ii) rotation alters the inverse energy transfer rate when Reλ≤255 and Ro≥0.8; (iii) the energy increase in large scale during inverse energy transfer exhibits a bell shape, the peak of which varies with Reλ and Ro